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paper, the expansion of [B. Gao, Phys. Rev. A 58, 4222 (1998)] is generalized to higher orders for 
angular momentum I > 2 to give a more complete description of ultracold scattering, including an 
analytic description of ultracold shape resonances of arbitrary /. We also introduce a generalized 
scattering length that is well defined and useful for all partial waves, to replace the traditional 
definition that fails for I > 2 due to the long-range van der Waals interaction. The results are 
used in a companion paper to derive analytic descriptions of atomic interaction around a magnetic 
Feshbach resonance of arbitrary angular momentum I. 
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I. INTRODUCTION 

One of the key factors that has made cold-atom physics 
such a thriving field has been the tunability of atomic in- 
teraction via Feshbach, especially magnetic Feshbach res- 
onances [U, 0, Q . While the s wave Feshbach resonances 
have received most of the attention for many years, the 
same concept is of course equally applicable to any other, 
nonzero partial waves. A number of such resonances have 
been observed experimentally 0, H, H, 0) 111 ; including 
a recent observation of a Feshbach resonance of / = 8 
Before we can understand questions such as what 
happens to a few-atom or a many-atom quantum system 
when a p wave or a d wave Feshbach resonance is tuned 
around the threshold, we need first to understand the 
corresponding two-atom system, namely, how to describe 
atom-atom scattering and atom-atom bound state with 
a Feshbach resonance of angular momentum I around the 
threshold. This is the subject of this and a companion 
study. 

Recall that for the s wave, a magnetic Feshbach reso- 
nance can be conveniently described by 



B - Bai-_ 



(1) 



Here ai^Q represents the s wave scattering length, which 
is tunable by the magnetic field, B, around a Fesh- 
bach resonance, a^gi^o is a background scattering length, 
Abi=o is a measure of the width of the resonance, and 
Bqi^q is the magnetic field at which ai=o — oo, corre- 
sponding to having a quasibound s state right at the 
threshold. Once the scattering length is determined, the 
scattering properties above the threshold, or the binding 
energy of an atom pair with large and positive scatter- 
ing lengths, can be determined, at least to a degree 
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from the effective-range theory (ERT) [nl, [T^, 13 1 which, 
for positive energies, corresponds to the expansion 



cot Si ^ -— + ^reik^ 



1 

ai 



(2) 
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where r^i is called the effective range. 

The situation for / ^ is considerably more complex 
due to the long-range van der Waals interaction, — Cg/r^, 
between atoms. First, we have the obvious problem that 
the scattering length is not defined for Z > 2 [ij, [lH . 
meaning that the energy dependence of the scattering 
phase shift will necessarily differ from that implied by 
Eq. (21). Second, a Feshbach resonance oi I ^ Q, when 
slightly above the threshold, is in fact a Fcshbach/shape 
resonance, the atoms in such a state see the angular- 
momentum barrier just like they would in a single chan- 
nel shape resonance state. When such a Feshbach res- 
onance is tuned from above to below the threshold, its 
width, due to the presence of the angular momentum bar- 
rier, can be expected to become increasingly narrow and 
approaches zero when it crosses the threshold to become 
a true bound state. Any complete theory for a Fesh- 
bach resonance oi I ^ Q has to be able to describe this 
increasingly rapid energy variation as it approaches the 
threshold. 

It is thus not surprising that the key to understand- 
ing Feshbach resonances of / ^ 0, to be presented in 
a companion paper, turns out to be the understand- 
ing of single channel shape resonances at ultracold en- 
ergies, or equivalently, the threshold behaviors of sin- 
gle channel shape resonances, which is addressed in de- 
tail in this work. We point out that this is a nontriv- 
ial problem that, to the best of our knowledge, has not 
been done except for the p wave description in Ref. . 
Part of the difficulty can be attributed to the breakdown 
of the semiclassical approximation around the threshold 

[n 113, in, Hi, no, nil, m m m, HI, which impues, in 

particular, that the tunneling amplitude or probability 
through the angular momentum barrier 25] cannot be 
obtained semiclassically. 
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Our analytic description of ultracold shape resonances 
is based on the small energy expansion of the quantum- 
defect theory (QDT) for — l/r^ type of interactions 
[3 HE HI]- It is presented here in a general analytic 
framework for ultracold atomic interaction that we call 
the QDT expansion, which is applicable with or without 
the presence of ultracold shape resonances. The presen- 
tation is organized as follows. In Sec. |TT1 we generalize 
the single-channel QDT expansion of Ref. [l^ to higher 
orders for I > 2 to give us a more complete description 
of atom-atom scattering around the threshold, includ- 
ing an analytic description of the threshold behaviors of 
single-channel shape resonances. In Sec. IIIIl we present 
further understanding of ultracold shape resonances by 
extracting from the QDT expansion analytic formulas 
for their position, width, and background. In Sec. lIVl we 
derive from the QDT expansion a generalized effective- 
range expansion, from which we introduce the concepts 
of a generalized scattering length and a generalized effec- 
tive range that are well defined for all angular momenta 
I. In Sec.|Vl we summarize the QDT expansions, derived 
previously in Ref. [131 > for the binding energies of the 
least-bound molecular state of arbitrary angular momen- 
tum using the standardized notations of Ref. ,25.] that 
we adopt here. We also present in this section a few inter- 
mediate results that will be useful in studying magnetic 
Feshbach resonances of arbitrary I. Further comments 
and remarks on the theory are presented in Sec. lVIl with 
conclusions given in Sec. lVIll The appendix presents an- 
alytic results of generalized scattering lengths and other 
QDT parameters, for arbitrary /, for two types of model 
potentials, a hard sphere with an attract tail of the type 
oi —Ce/r^, and the Lennard- Jones potential of the type 
of LJ(6,10). 

The analytic description of atomic interaction around a 
magnetic Feshbach resonance of arbitrary I , which is nec- 
essarily multichannel in nature P, |3i @] j is developed in a 
companion paper. It will be accomplished through first, a 
rigorous reduction of the underlying multichannel prob- 
lem, as described by the multichannel quantum-defect 
theory (MQDT) for type of interactions [28], to an 

effective single-channel problem, and a subsequent appli- 
cation of the results of this work. 



II. QDT EXPANSION FOR SINGLE-CHANNEL 
SCATTERING 

In this section, we derive and discuss the QDT expan- 
sion for single-channel scattering. It differs considerably 
from the ERT which assumes that fc^'+i cot Si 

is both an analytic and a slowly varying function of en- 
ergy around the threshold, and can therefore can be ap- 
proximated by the first few terms in its energy expan- 
sion. The QDT expansion makes no such assumptions. 
The only quantities expanded are the universal functions 
of a scaled energy that are associated with the long- 
range potential ^25j. There is no assumption about how 



fc^'+^cotJ; or tan (5; may vary with energy. It is partly 
for this reason that the QDT expansion can give ana- 
lytic description of an ultracold shape resonance, which 
has rapid energy dependence that would have required, 
at least, partial summation over all orders of energy in 
more standard approaches. The QDT expansion is also 
more than a small-energy expansion. It is simultaneously 
a large-/ expansion, in the sense that for any fixed energy, 
there is a sufficiently large / beyond which it becomes ap- 
plicable. As a related consequence, the energy range over 
which the QDT expansion is applicable increases rapidly 
for larger I. 

While the derivation of the QDT expansion is some- 
what tedious, the end result will be very simple. It is 
given by a single analytic formula applicable to all angu- 
lar momentum / and regardless of whether there is, or is 
not, an ultracold shape resonance. To make the deriva- 
tion easier to follow, we first rewrite, in Sec. Ill Al the 
QDT for single channel scattering [l^, [H, [2^ in a form 
that makes its subsequent expansion, presented and dis- 
cussed in Sec. Ill Bl fully transparent. 

A. QDT for single-channel scattering 

For a single channel problem with long-range —Ce/r^ 
interaction, the scattering above the threshold is de- 
scribed rigorously in the QDT by the following equation 
for the K matrix [H, HI ill] 

Ki = tan 5, = {Z^^.K^ - Z%){Z% - Zl^K^y' . (3) 

Here K'^{e,l) is a short-range K matrix that depends 
weakly on both the energy e and the angular momentum I 
[2^. The Z'^y are elements of the Z"^ matrix for the — 
type of potentials, labeled here using the standardized 
notation of Ref. [25] . They are given explicitly by 

Z% = Af { [1 ~ {-ifM.j tan^(z/ - v^)] sin(7ri./2)X, 
+ [1 + (-l)'Me,ztan^(z/- z/o)] cos(7ri./2)yz}(,4) 

Z% = Af { [tan7r(i. - uq) - (-I)'M,,,] sin(^z//2)Xz 

+ - z/o) + (-l)'M,^,] cos(7ri//2)rz} , (5) 

Zl^ = Af { [1 + (-l)'M,^z tan^(i/ - v^)] cos{TTiy/2)Xi 

- [1 - (-l)'Me,ztan^(z/- z/o)] sm{7n^/2)Yi}(,6) 
Zl^ - Af { [tan7r(i/ - v^) + (-I)'M,^,] cos(7ri//2)X, 

- [tan7r(i. - v^) - (-l)'A/e,;] sin(W2)y,} , (7) 

where 

j^z _ GeA^) cos7r(t/ - vq) 
' " \/2(X2 +y,2)sin7ri/ ' 

1^0 = (2/-t-l)/4, Me^i = Ci(-z/)/G'e^z(i'), with the char- 
acteristic exponent and functions Xi, YJ, and G^^i be- 
ing defined in Ref. [29|. They are all universal functions 
of a scaled energy, 

es = e/s_E, (8) 
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TABLE I: Sample scale parameters for A + A type of systems 
where A is an alkali-metal atom. The /Sg — (2^C(,/'h?Y^'^ is 
the length scale. The se ~ {Tl^ /2fi){l/ (3^)^ is the correspond- 
ing energy scale. It is given both in units of /uK and in units 
of MHz. St ~ h/sE is the corresponding time scale. All are 
determined by the Ce coefficient and atomic masses. 



Atom Ce (a.u.) 


(a.u.) 


SE/ks {^J.K) SE/h (MHz) 


St (ns) 


''Li 1393.39" 


62.52 


7368 


153.5 


1.037 




89.86 


933.1 


19.44 


8.186 


*°K 3897' 


129.8 


257.3 


5.360 


29.69 


**^Rb 4707" 


164.3 


75.58 


1.575 


101.1 


^^^Cs 6860"^ 


201.9 


31.97 


0.6662 


238.9 




"From Ref. [30 
''From Ref. i3l|| 
'^From Refs."0,l33 
"^From Ref. E 



where se = (?i^/2/^)(l//36)^ is the energy scale corre- 
sponding to the length scale Pq = {2fiCQ/h'^y^^ that is 
associated with the —C^/r^ potential. For the purpose 
of providing orders of magnitudes, these scales, and also 
the related time scale, st = 'hj SE^ are tabulated in Ta- 
ble [T] for selected alkali-metal dimers. They are all de- 
termined by the Cg coefhcient and the reduced mass [i. 
Equation ([3]) gives an exact description of single-channel 
atomic scattering in terms of a set of universal functions 
determined solely by the long-range interaction. All the 
short-range physics are encapsulated in K'^{e, I), which is 
slowly varying in both e and I around the threshold [1^ . 

Instead of the K'^ parameter, the short-range atomic 
interaction in any partial wave can also be described us- 
ing alternative parameters such as the quantum-defect 
/i'^ or the parameter Kf^ [^B]. They are all related, but 
have different utilities for different purposes. For this 
work, which focuses on the ultracold regime around the 
threshold, the most convenient parameter is -flTf^. It is 
defined as the short-range K matrix associated with the 
and g'^^ reference pair of Ref. [2^, and is related, for 
— type of potentials, to the K'^ and fi'^ by 



K^ieA) - tan(7riyo/2) 
1 -htan(7ri^o/2)-ft:'=(e,0 
tan[iTif{e, I) - In /A] . 



(9) 
(10) 



With this definition, if; (e = 0) = corresponds to hav- 
ing a bound or quasibound state of angular momentum I 



right at the threshold, a small and positive Ki^{Q) corre- 
sponds to having a bound state of I close to the threshold, 
and a small and negative Kf^ifi) corresponds to having 
a shape resonance of angular momentum I close to the 
threshold [23, jH, [3^ 
called xi in Refs. [27L [3 

Ref. [H by Kf = -Kf, allows it to be used as an expan- 
sion parameter when there is a state close to the threshold 
[131. As will become clear later in the paper, is also 
the short-range parameter that has the simplest relation 
to the scattering lengths and the generalized scattering 
lengths. 



This property of Kf , which was 
and is related to the K^* of 



Defining 9i by 



sm Ui 



cos Bi = Xi{Xf + Yi 



-1/2 

2^-l/2 



(11) 

(12) 



the elements of the Z'^ matrix can be further written as 



sm [ — TTt/ + f; 



-{—l) M^j tan7r(iy — Vq) sin —t:i> — 9i 



tan7r(i^ — uq) sin ( —TTiy + 9i 



cos I -^TXV + d[ 



-|-(— 1) Afe,; tan7r(i^ — Vq) cos [ -tti^ — Bi 



where 



tan7r(i/ — vq) cos [ 2""*^ ^" ^' 



+ (-l)'Me,,C0S ( -TTV-Bi 



,z _ Ge^;(i/) cos 7r(zy - 



Bf = 



(;3) 



(14) 



(15) 



(16) 



(17) 



Substituting Eqs. (IT3l)-([T6l) into Eq. ([3]) and using Kf{e) 
as the short-range parameter, we can rewrite the QDT 
equation for the K matrix as 



J 



Ki = tan (5; = - tan[7r(zy - z^o)] - Ais{es)k 



21+1 



1 — tan^[7r(i/ — j/q)] 



1 + Ai,{es)kf+' tan[7r(j/ - 1^0)] 



(18) 



where = e 



1/2 



fc/?6, and 



Als{€s) 



[1 + {-iyKf]{l + tan 61, tan[7r(i/ - i/o)/2] - [(-1)' - Kf]{i&n6i - tan[7r(j/ - J/o)/2]} 
Kf - tan 6*, - tan[7r(z/ - i^o) /2] - Kf tan 6*, tan[7r(i/ - v^) /2] ' 
I 09) 



This is stih an exact expression for Ki, and it has been 
written in a way to make its QDT expansion fuUy trans- 
parent. 



B. QDT expansion 



rive straightforwardly from the analytic solution for the 
-l/r^ type of potential [29t : 



6^ + 0(4), (20) 



For small energies around the threshold, or for arbi- 
trary energy but sufficiently large the quantities in 
Eq. (fT8|) can be represented by expansions [23| that de- 



tan 6] = — 



22(^1-1) 



(21) 



24-2-1 sin(7r!.o) [r(i.o)r(i.2 + 1)]^ 



X h +2(j/-z.o)ln|e.| 



1 



26 



1 



(-l)'7r + 81n2 + — + 27/'(i/o) 
16 1 



+0[|e.|-+4(ln|e.|)2] 



r 



(22) 



where vi^ = (2Z + l)/4, as defined earlier, V2 = 2vq = Combining these expansions with Eqs. (US]) and ([TO]) 

Z -I- 1/2, and ip{x) is the digamma function [s^l- shows that the X matrix has the following structure 



tan (5, — —i:{v — vq) — asik 



21+1 1 



Kf 



01 - tt{j^ - j^o)/2 + 0(fcf ) 



[1 + 0{kt In ks)] + 0{k!) , 

m 



where 



asi 



24'+i[r(//2-K l/4)r(? + 3/2)]2 

TT 

22'-i[r(V2 + l/4)]2[(2Z-Hl)!!]2 ' 



(24) 



is what we call the scaled mean scattering length for an- 
gular momentum /. It is a generalization of the s wave 
mean scattering length of Gribakin and Flambaum [l9| 
to p [131 and higher partial waves. More explicitly, 



2tt 



asi=o — 



asi=i 



[r(i/4)]2 
[r(i/4)]2 



367r 



0.4779888 , 



0.1162277 . 



(25) 
(26) 



For larger /, it decreases rapidly and can be either com- 
puted directly or obtained from asi=o and asi=i using the 
recurrence relation 



asi+2 



[{21 + 5)i2l + 5){2l + 1)] 



(27) 



which gives, e.g., asi=2 = asi=a/225. Together, they rep- 
resent a set of universal numbers determined solely by the 
long-range potential. The corresponding mean scattering 
length, with scale included, is defined as di = asi(3l^'^^ . 

Equation ((23)) . if kept to the order of fcf, gives the 
following result for the s wave 



5 



tan 



1=0 



15 



Kf 



{l-Kf){lkl + f^ki) 



Kf 



15 



:fc„ In ks 



2 

15 



22 



In 2 - 7 fc; 



(28) 



where 7 = 0.5772156649 ... is the Euler's constant [37|. 
Resuhs of even higher orders are possible. An s wave 
result to the order of fcf was given in Ref. [lB| • Similar 
results can be written down for other partial waves. 

For the range of energies of interest in cold-atom 
physics, this level of complexity is actually unnecessary, 
especially for higher partial waves. For simplicity and 
easier application, we will use the following approxima- 
tion of Eq. ^ for aU I 



tan (5; « K, 



(B) 



(D) 



(29) 



TABLE II: Critical scaled energies, esd, for different angular 
momentum I. Beyond esd, the characteristic exponent f for 
the — solutions [2^ moves off the real axis and becomes 
complex. 



I 



where 



1.544707 
2.358067 
6.891073 
25.29322 
33.17273 
89.85261 



6 
7 
8 
9 
10 
11 



94.82401 
212.3067 
206.8228 
406.8301 
383.5433 
687.4041 



{B) 



-7r(j/ - vq) 



{21 + 5){2l + 3){2l + 1){21 - 1){21 - 3) 



£^(30) 



is a term that could have been derived from the Born 
approximation (see, e.g., Ref. [Ill), and 



(D) 



21+1 



(31) 



describes the deviation from the Born term. Here 



asi 



(-1)'+ 



1 + KfOi 



Kf 



{2l + i){2l-l)- 



Kfe., 



{21 + S){21 - l)Kf + e, + wie 



(32) 



3) 



can be regarded as a scaled, energy-dependent, general- 
ized scattering length. Both expressions for Asi{es) are 
useful for different purposes. In Eq. ((33|) . wi is an l- 
dependent constant defined by 



Stt 



Wl 



2{2l + 5){2l + 1){21 ~ 3) 



(34) 



Equation p2|) . while less exphcit compared to Eq. ((33)) . 
is more convenient for a number of conceptual purposes. 
Here 7r(i' — i^o) is given by Eq. (ISO]) , and 9i is given by 
Eq. (I2ip . or more explicitly by 



9i 



1 



(2/-f 3)(2/- 1) 



(35) 



Together, Eqs. give a single analytic formula 

for low-energy atomic scattering that we call the QDT 
expansion. It is applicable to all with or without the 
presence of ultracold shape resonances. It has the follow- 
ing additional characteristics, (a) The QDT expansion 



remains applicable no matter how rapidly the Kf{e) pa- 



rameter may depend on energy. In arriving at Eqs. 
([33]) . the only quantities that are expanded are the uni- 
versal functions associated with the long-range potential. 
There is no assumption about the values of iff°(e), or 
how it may depend on energy. While this feature of the 
QDT expansion is not important for true single channel 
cases, for which the energy dependence of Kf is almost 
always negligible [H, HI] , it will become crucial when we 
apply it to effective single channel problems that are de- 
rived from intrinsically multichannel ones, for which the 
energy dependence of Kf is generally important, (b) 
The energy range over which the QDT expansion is ap- 
plicable increase rapidly with /, roughly as l^. This is 
consistent with the fact that the QDT expansion is si- 
multaneously a large I expansion [27[. For any energy, 
there is a sufficiently large I beyond which it becomes 
applicable. The only condition for the applicability of 
the QDT expansion is <C esd, where Csd is the scaled 
critical energy beyond which the characteristic exponent 
becomes complex [29^ . Values of e <;c/ for the first twelve 
partial waves are listed in Table [Til For large I, they cor- 
respond roughly to the scaled height of the angular mo- 
mentum barrier, given for — 1/r^ type of potentials by 
Hd = (2/33/2)[;(/ -I- l)]3/2. As win be discussed in more 
detail elsewhere, the critical energy is also the energy 
around which the behavior of the system goes from being 
quantum to being semiclassical [23, 25]. Thus the QDT 
expansion can also be regarded as a quantum expansion, 
applicable over a quantum region of energies where the 
quantum reflection probability is close to 1 [29| . 

Figures[l]and[2]illustrate the QDT expansion for p and 
d waves, respectively, using cases that have a shape reso- 
nance in the threshold region. Since our interest here is in 
the case of a single channel, Kf is taken to be a constant 



\m,m 



Kf 



Kf{e = 0), which is related in a simple 
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FIG. 1: (Color online) p wave partial cross sections for 
ai=i/ai=i = -20 (corresponding to iCfi^i « -0.0526). Re- 
sults of the QDT expansion (dashed line) are compared to 
the exact QDT results computed from Eq. (|3]) (solid line), 
and the results of ERT (dash-dot line). 



FIG. 3: (Color online) s wave partial cross sections for 
ai=o = 0, representing the worst-case scenario for the QDT 
expansion. The results of the recommended QDT expansion 
(dashed line) are compared to the exact QDT results from 
Eq. ([3} (solid line), and the results of a higher order expan- 
sion, given by Eq. (|28|l (dash-dot line). 
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FIG. 2: (Color online) d wave shape resonances as described 
by the QDT expansion for different values of Ki2,2- The ex- 
act QDT results computed from Eq. ([Sj are only shown for 
is:f°2 = -0.04 as they become indistinguishable from those of 
QDT expansion for smaller li(';'i22l- F™' higher partial waves, 
the QDT expansions are more accurate, and applicable over 
a greater range of scaled energies than shown here for the d 
wave. 



way to the scattering length by Eq. (|l8|) or ([521) . to be 
discussed in more details later. Note that the QDT ex- 
pansion is applicable regardless of how narrow the shape 
resonances may be. In fact it is more accurate for nar- 
rower resonances which are necessarily located at smaller 
energies, as illustrated in Fig.O For higher partial waves, 
the results of the QDT expansion are no longer distin- 
guishable from the exact QDT results, computed using 
Eq. ([3]), in the range of energies shown in Figs. [T] and 
[21 Such results are easily calculated from the analytic 
formula, and are therefore not shown. 



Figure [Tjis also used to illustrate the limitation of the 
ERT nil, 111 [11. In the figure, the ERT results are 
computed from [see Eq. ^] 



cot 5i= 



-l/ai-_ 



(36) 



Note that since the effective range is not defined for the 
p wave due to the van der Waals interaction [TJ, Ul , it is 
incorrect to add an effective range term to Eq. ([55[) . The 
lowest order correction to the right-hand-side of Eq. p6[) 
is of the order k (see Eqs. ([in|)-([23|), or Ref. [H), not of 
the order of k^ as implied by ERT. Figure [T] clear illus- 
trates the failure of ERT, which misses the p wave shape 
resonance completely. The consequence of such a failure 
for two atoms in a trap has been illustrated elsewhere 
[3^ . Careful readers should note that ERT description 
of p wave interaction is still sometimes incorrectly used 
in the literature. This is unfortunate, considering that a 
correct description, at least in the case of a single chan- 
nel, has been available for some time |T5| . 

Among all partial waves, the QDT expansion for the s 
wave has the smallest range of applicability because the 
critical scaled energy, given in Table [Til is smallest for 
I — 0. For the s wave, the most stringent test of the 
expansion is for ai^o — 0. It corresponds to a case where 
the lowest order contribution, —ai^^k, to tan 5; goes to 
zero, and we are directly testing the higher order terms. 
As will be discussed in more detail in Sec. lIVi this is also 
the case where the effective range expansion, even in its 
generalized form of Sec. llVi fails completely. Figure [31 
illustrates the accuracy of the QDT expansion for this 
worst case scenario. It shows that the simpler QDT ex- 
pansion that we are recommending here is accurate for 
roughly eg < 0.1. From Table [11 it is clear that this range 
of scaled energies already covers all energies of interest 
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in cold-atom physics. The QDT expansion is more accu- 
rate and applicable over a greater range of energies in all 
other cases and for all other partial waves. 

For the s and p partial waves, our results here are 
consistent with those of Ref. except they are now 
written in simpler forms that are also better for physical 
interpretation. For / > 2, the if term, which is of 
the order of fc^'+^ under nonresonant conditions of |A'f"| 
being of the order of 1 or greater, is normally negligible, 
as was done in Ref. llSl. In doing so, however, we missed 
in our previous work [15| the analytic description of low- 
energy shape resonances for Z > 2, and an opportunity to 
define the generalized scattering length and the general- 
ized effective range. These subjects are addressed in the 
next two sections, respectively. 




III. THRESHOLD BEHAVIOR OF SHAPE 
RESONANCES 

The QDT expansion of Eqs. (HHl-dSSl) is applicable 
whether or not there is a shape resonance in the thresh- 
old region. When such a resonance does exist, which 
occurs for I > 1 and a small and negative K^'^, namely 
for iff^ < and |iff°| ^ 1, further conceptual under- 
standing can be achieved by extracting from the QDT 
expansion the standard parameters characterizing a res- 
onance, namely its position, width, and background (see, 
e.g., Ref. [401). 

The position of a shape resonance in the threshold re- 
gion can be determined from the root of the denominator 
in Eq. (|3T|). 



FIG. 4: (Color online) Universal relation between the scaled 
lifetime and the scaled position of a single-channel shape reso- 
nance, as given by Eqs. (|4ip and (|45p . for — 1/r'' type of long- 
range potentials. The lifetime can change by many orders of 
magnitude for a small change in the resonance position. 



^bgi'i'^s), given by 



bgl = l ~ 1 _ iLc , 



' 3/2 3/2' 
— ^sl 



14 



21a, 



'5/2 5/2- 
^sl 



1=2 



1 + J^^sl \ (-S — ^sl 



1 1 a,,63/2 , 
(42) 

(43) 



9i + Tr{iy - i/o)/2 



(37) 



It has solution only for a small and negative iff", for 
which it can be solve perturbatively to give the scaled 
resonance position as 

esi w -(2/ + 3)(2; - l)Xf° [1 -I- wi{2l + 3)(2Z - l)Xf°] , 

(38) 

in which wi is defined earlier by Eq. p4p . Around the 

shape resonance, the term if/^'' = —Asi{es)k'f^^ can be 
written as 



^bgl y'^s) 



2 Cs — £«; 

with the scaled width, 'jsi, given by 

jsi - 2[(2? + 3)(2/-l)]'+3/2a,;(_iff")'+i/2 
X 1 



(39) 



3^(2? + 3)(2;-l)^,o 



4(2/ + l)(2Z-3) ^ ' 
2(2/ + 3)(2/-l)a,;(£,;)'+i/^ 
1 -I- 2wiesi 

(D) 



(40) 
(41) 



and the contribution of Kl to the background 



{D) 



l>2 







(44) 



Corresponding to the width of Eqs. (j40|) and (|4T|) . there 
is a well-defined hfctime (see, e.g., Ref. [40l |) 



n = TslST = il/jsl)sT 



(45) 



where St = T^/se is the time scale associated with the 
length scale /Jg, with sample values for alkali-metal atoms 
given in Table [H and Tsi = ti/st — ^/"fsi is the scaled 
lifetime. 

Equations and ([1T|) show that the width of a 

shape resonance goes to zero as it approaches the thresh- 
old, below which it becomes a true bound state. The 
corresponding lifetime goes to infinity. Equation (HT|) 
further shows that the scaled width, and therefore the 
scaled lifetime, when viewed as a function of the scaled 
resonance position, follows a universal behavior that is 
uniquely determined by the long-range interaction. In 
other words, while the resonance position itself, as given 
by Eq. ([55]) . depends on the short-range parameter 
the functional form of the scaled width or lifetime versus 
the scaled resonance position is independent of it. This 
universal behavior is illustrated in Figures S] for the first 
few partial waves. It also illustrates that the lifetime of 
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a shape resonance can change by many orders of magni- 
tudes, from a microscopic scale of st (see Table H]), to a 
macroscopic scale of seconds or longer, as the shape res- 
onance approaches the threshold. The existence of such 
a potentially macroscopic time scale ,9'] is one of the key 
differences between the coupling of atoms in ^ 7^ and 
I = partial waves. 



length for an arbitrary I. This, and the definition of a 
generalized effective range for an arbitrary I, can be done 
in a way that more closely resembles to the standard ERT 
[H [Hill, as follows. 



IV. GENERALIZED SCATTERING LENGTH 
AND GENERALIZED EFFECTIVE RANGE 



The QDT expansion, Eqs. (|29|) - ((33|) . already suggests 
that one may be able to define a generalized scattering 



Defining s'j-'^^ 



-^{v — i/q), we have from Eqs. ([18] 



Kf -61+ 0{ki 



asij 1 



[l + 0(fc^lnA:,)] +0{kt) 



(46) 



In the case of a single channel, the energy dependence of 
is completely negligible [l^, [2^. A further expan- 
sion of the denominator in Eq. (l46l) gives the following 
generalized effective range expansion 



cot(,5i - 5 



(B) 



-^ + lr,ik^ + 0{kHnk) , (47) 
ai 2 



with the generalized scattering length, a;, given by 



ai = ai 



1 



Kf{e = 0) 



(48) 



where a; = a^i/Jg'^^ is the mean scattering length for 
angular momentum / that we have defined earlier (with 
scale included). The generalized effective range, r^i, is 
given by 



rei 



{2l + 3){2l-l) 



^ {1 +[(-!)' -(a,M)]^} 



(49) 

For quantum systems with a long-range type of 

interaction, Eq. (l47l) defines the generalized scattering 
length and effective range for an arbitrary I. It coincides 
with standard definitions of scattering lengths and effec- 
tive ranges whenever they are well defined in the standard 
theory [ill, [H, El, El . Namely, 5/ = a/ for ^ = and 1, 
and Tel = Tel for ^ = 0. Similar to the standard theory, 
a; has a dimension of /3g'^^, and Vei has a dimension of 
p-2i+i^ The generalized scattering length is also related 
to Asiics) by ai = I,i(e, = 0)Pf+\ 

With this definition of generalized scattering length, 
having a bound or quasibound state right at the thresh- 
old, characterized in terms of the K^'^ parameter by 
^I'^i^ = 0) = 0, always corresponds to ai = 00, for any I. 



Similarly, having a shape resonance close to the thresh- 
old, characterized by a small and negative -fi'f*', corre- 
sponds to having a large and negative generalized scat- 
tering length; having a bound state close to the threshold, 
characterized by a small and positive iiTf", corresponds to 
having a large and positive generalized scattering length. 
Such similarities to the s-wave interaction make the gen- 
eralized scattering length an easy parameter to under- 
stand, without having to know the QDT behind it. As 
examples of the generalized scattering lengths for arbi- 
trary I, we present, in the Appendix „ their analytic re- 
sults for two classes of model potentials with — 1/r^ type 
of asymptotic behaviors. 

We point out that the main utilities of the generalized 
effective range expansion, Eq. (|47)) . are (a) to define the 
generalized scattering length as an alternative parame- 
ter for describing low-energy atomic interactions, (b) to 
make a connection between the QDT expansion and the 
ERT to the degree possible, and (c) to simplify the un- 
derstanding of the QDT expansion for peoples who are 
not completely comfortable with QDT formulations. It 
is not meant to be a replacement for the QDT expansion. 
As far as accuracy is concerned, the QDT expansion is 
always more accurate. This loss of accuracy in the gen- 
eralized ERT occurred in expanding the denominator of 
Eq. (j46p . whose full representation would have required 
an infinite number of terms in the standard ERT type of 
expansions. 

The procedure of expanding the denominator has more 
severe consequences in the special case of a; = 0, for 
which Vei = 00 from Eq. (|49|) . and the effective range 
expansion, even in its generalized form here, becomes 
meaningless. In comparison, the QDT expansion remains 
applicable, and gives, for a; = [corresponding to Kf'^ = 
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-i-iy 



K 



{2i + 3){2i-i)-{-iyk^-{-iywikt 



(50) 



This result also implies that a; = changes the threshold 
behavior of Ki = tan (5/ for ^ < 2. For a; 7^ (and a; 7^ 
00), it is clear from the QDT expansion that the threshold 
behavior for ^ < 2 is determined by the K^'^^ term that 

behaves as K^'^^ ^ —aik^'-^^. The threshold behavior 
for ? > 2 is dominated by the Born term, which behaves 
as kg. For a; = 0, Eq. ((50)l means that the threshold 
behavior for the s wave changes from tanSi^Q ~ —ai=ok 
to tan(5;=o ~ —'^o,si=ok^ /3. The threshold behavior for 
the p wave changes from tan(5;=i ^ — aj^ifc^ to being 
dominated by the Born term tanSi^i — 7rfc^/35. 

Another special case for which the threshold behavior 
may be modified is the case of a; —00 {K['^ — 0), cor- 
responding to having a bound or quasibound state right 
at threshold. The QDT expansion gives for this special 



(2/ + 3)(2r 



Da.kf-' 



1 + W/fc^ 



It means, in particular, that the threshold behavior for 
tan Si is modified for I < 2 for having a bound or qua- 
sibound state right at the threshold. Specifically, it is 
changed from tan (5/ ~ — aifc^'+^, for cases of a; 7^ 00 (and 

ai ^ 0), to tan (5; {21 + 3){2l - l)dsikf-^ for a; = 00. 

The threshold behaviors for ^ > 2 remain dominated by 
the Born term (~ k^) even with a bound state right at the 
threshold. We note that the generalized effective range 
expansion would have given = — 2/3|/[(2/-|-3)(2Z— l)a/] 



and tani5; 



r(B) 



{21 + 3)(2Z - l)asikf-^ for ai = 00, 



corresponding to the lowest order term in Eq. (|5ip . 

As a further comment on the generalized effective- 
range expansion, we note that in arriving at Eq. (|49p for 
the effective range, we have assumed that the energy de- 
pendence of is negligible. This means that Eq. (HS)) 
is, strictly speaking, a single channel result that will need 
to be modified for effective single channel problems where 
the energy dependence of Kf'^ is generally important (to 
be discussed in detail in the companion paper). With 
this limitation in mind, Eq. (|49p does imply that in the 
case of a single channel, for which the energy dependence 
of is negligible, the generalized effective range is not 
an independent parameter, but can be determined from 
the generalized scattering length. It further implies that 
the generalized effective range has the property of r^i > 
for I — and rei < for I > 1, again rigorous only for 
true single channel problems. 

All QDT expansion results of previous two subsections, 
which were parametrized using i^f°, can be written in 



(52) 



terms of a; using Eq. (j48p . or equivalently. 



(e = 0) 



1 



We do not give these expressions explicitly, in part to 
again emphasize the following subtle, but important 
point. The expressions in terms of iff*^ are more gener- 
ally applicable because they make no assumption about 
the energy dependence of K^^. The corresponding ex- 
pressions in terms of a/ automatically assumes the weak 
energy dependence of i^f°, since in effect we are using 
Ki^{e = 0) at other energies. This subtlety is not an 
issue for true single-channel problems, but becomes one 
for effective single channel problems derived from multi- 
channel cases (see the companion paper), for which the 
equations in terms of Ki'^ remain applicable, but gener- 
ally not those in terms of a; . 



V. QDT EXPANSION FOR BOUND STATE 
ENERGY 

The bound spectrum of a two-body, single channel sys- 
tem with — l/r^ type of long-range interaction is given 
rigorously by the solutions of 1^, 2^, [2^ 



where 
X?(e.) 



xU^s) ^ K^{e,l) . (53) 

tan6'i tan(7rj//2)(l + M^j)/{1 - M,j) 
1 - tan6'; tan(7rjy/2)(l + M^j)/{1 - M^j) 



(54) 

is a universal function of Cs that depends only on the 
exponent of the long-range interaction and on the angular 
momentum I. 

For deriving the QDT expansion for the energy of the 
least-bound state that is close to the threshold, it is again 
more convenient to rewrite Eq. (|53p in terms of the Kf^ 
parameter, as 



where 



x!'{^s) 



xf{es)^K^"{e), 

_ xF(es) - tan(7r;/o/2) 
~ l + tan(7ri/o/2)x,(es) 



(55) 



(56) 



Using expansions of v, 9i, and M^j, as given by Eqs. ((20|) - 
(P^ . we have, for sufficiently small energies or sufficiently 
large I, 



% + Ti{i' - ;/o)/2 -I- asi=oKs - asi=oKsdi 
1 — asi^QKs — asi=oK,s6i 



1 - asi=iKl 



x'i>2{^s)~ei+n{u-v^)/2 



(57) 



(58) 



(59) 



Here = (— e.,)^/^, and 9i and 7r(i^ — z^o) £^re given by 
Eqs. ([55)1 and ([50]) . respectively. 
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Having a bound state of angular momentum I close to 
the threshold corresponds to having a small and posi- 
tive Ki ^. The energy of such a state can be obtained by 
solving Eq. (j55p perturbatively using expansions given 
by Eqs. j^-dM]). This has been done in Ref. We 
summarize the results here using the standardized no- 
tation {xi of Ref. prj is renamed Kf°), for the sake of 
completeness and easy reference. 

For the s wave, we obtained 



energy variations of i^f", will be useful in developing the 
analytic description of magnetic Feshbach resonance of 
arbitrary I, to be presented in a companion paper. 



VI. DISCUSSIONS 



esi=o — — ^ 



1=0 



We make here some miscellaneous comments on vari- 
ous aspects of the theory, which, while not essential, may 

(iirf'^g)^ [l -I- giKi^Q + g2iKi^Q)^]+0{{K[^Q)^ helpful in the understanding of this and related theo- 
~ _ _ _ ^^^^^^ 

(60) (a) The QDT expansion breaks naturally into two 



where e^; is the I wave bound state energy scaled accord- 
ing to Eq. dH]), and gi = 2[l/(3a2;^o) - 1] « 0.9179195, 
and c/2 = (5/4)5? - 2 w -0.9467798. 
For the p wave, we obtained 



cO 



1 



xl/2 



cO 
1=1 



1.299466, and ho 



(61) 
= 2,hl/2 - 



where hi = 5'^/^as;= 
57r/14« 1.410919. 

For / > 2, the energy of the least-bound molecular state 
is given by Eq. ([55]) . namely the same equation that gives 
the position of the shape resonance above the threshold. 
This comes from the fact that Eq. ([55)1 , with xf' given by 
Eq. ([Sg]) . is the same as Eq. ([37|) . which determines the 
shape resonance positions. The only difference is that 
while Kf' is small and negative for a shape resonance 
close to the threshold, it is small and positive for a bound 
state close to the threshold. 

In cases of true single channel problems, these equa- 
tions can again be written in terms of the generalized 
scattering length using Eq. ([5^ . The corresponding 
equations for the s and p waves can be found in Ref. [13] . 
We skip writing down these equations explicitly to again 
emphasize the greater applicability of the equations in 
terms of . Namely, unlike the equations in terms 
of the scattering length or generalized scattering length, 
which assume the energy-independence of Kf^, the equa- 
are, strictly speaking, applicable 
itself depends on energy. The only sub- 
tlety here is that when the right-hand-sides of equations 
such as Eqs. (|60)) and ([6T|) become energy-dependent, 
they need to be solved again to obtain the bound state 
energies. 

The results presented in this section have been verified 
in Ref. [23| through comparisons with the exact QDT re- 
sults, obtained by solving Eq. ([55)1 or (|55|) numerically. 
The result for the s wave, Elq. ([50]) . has also been veri- 
fied by Derevianko et al. |41| and by Julienne and Chin 
p3 | through independent numerical calculations. In sum- 
mary, the formula for the s wave bound state energy is ap- 
plicable for ai=o > 2a;=0i corresponding, approximately, 
to ai^o > Pa or |es;=o| < 4. The formulas for other par- 
tial waves are applicable over a greater range of bound 
state energies |27|. The intermediate equations, including 
the expansions of the function as given by Eqs. ([F7)) - 
and Eq. which is applicable regardless of the 



tions in terms of K^^^ 
even when 



terms, a Born term , and a deviation from the 

Born term k[^\ This separation of Ki = tan 5/ is also 
convenient for a number of other applications including, 
e.g, the understanding of angular distribution (see, e.g., 
,Refs. (43ll43|V There are other separations of Ki ~ tan 5/ 
possible. For example, it can be separated into a back- 
ground term and an interference term, as suggested in 
Ref. [1^. One can show, however, that these two types 
of separations are the same for / > 2, as is reflected in 
the fact that the K^'^'' term has negligible contribution 
to the background for / > 2 [see Eqs. (|39l) and (l44l) ]. 

(b) Unlike the cases of s and p partial waves, the gen- 
eralized scattering length for ^ > 2 has little effect on 
scattering cross sections around the threshold, except in 
determining the location of the shape resonance, and if 
there is one, the cross sections around it. 

(c) For true single channel problems, having an ultra- 
cold shape resonance can only happen by accident. In 
this sense, the most important applications of a theory 
for ultracold shape resonance are not to the true sin- 
gle channel problems, but to the effective single channel 
problems associated with Feshbach resonances, for which 
there will be at least one state in the threshold region as 
the Feshbach resonance is tuned around it. It is for this 
reason that we kept stressing the formulas that remain 
applicable when Kf*^ becomes energy-dependent, which 
is the single most important difference between an ef- 
fective single channel problem for a Feshbach resonance, 
and a true single channel problem. 

It is for the same reason that we did not emphasize the 
relationship of JsTf" for different /. For a true single chan- 
nel problem, the Kf^ parameter for different / are not in- 
dependent, but are related to each other through Eq. 
in which K'^ is approximately independent of / (26,] . As 
a result, the Kj^^ for the first few partial wave can all be 
determined from, e.g., the s wave scattering length. This 
simple relationship between different I breaks down for 
the effective single channel problems for Feshbach res- 
onances. For such intrinsically multichannel problems, 
the Kf'^ parameters for different / are still related, as 
the underlying short-range K'^ matrix is still approxi- 
mately /-independent. Their values for different I can 
still be computed from, e.g., the singlet and the triplet s 
wave scattering length for alkali-metal atoms [1^. Their 
relationship, however, becomes more complicated than 
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Eq. ([9]) and less transparent. 

(d) While we have no intention of promoting one pa- 
rameter over the other, as they all have different utilities, 
we hope this work again illustrates that the scattering 
length, or the generalized scattering length, is not the 
only parameter, or necessarily the best parameter for 
characterizing atomic interaction at low temperatures. 
While many of our results, such as those for resonance 
positions and binding energies, can be written in terms of 
the scattering length or the generalized scattering length, 
such equations are more complex than those in terms of 
iff", and have more restricted applicability to the case 
of single channel. Furthermore, it should be clear that in 
the scattering length representation, what is important 
is not the scattering length itself, but the dimensionless 
ratio ai jai , as is also recognized by Chin et. al Q for the 
s wave. 

(e) We point out that there are a number of other 
successful QDT formulations of atomic interaction that 
are conceptually similar to ours in maiiy aspects, but 
based on numerical reference functions [3, |4^, [Ifl 
[isl . |49| . Without incorporating at least some ingredients 
of an analytic solution, numerical solutions generally run 
into difficulty for very large scattering lengths, namely 
when there is a state very close to the threshold. This is 
true for the s wave [50| . and more so for higher partial 
waves. 



VII. CONCLUSIONS 

In conclusion, we have presented an analytic descrip- 
tion of atomic interactions at ultracold temperatures for 
the case of single channel. In particular, a QDT expan- 
sion for scattering has been developed that is applicable 
to all angular momentum with or without the presence 
of ultracold shape resonances. Using the QDT expan- 
sion, we have developed a fully analytic characterization 
of ultracold shape resonances in terms of its position, 
width, and background. We have also introduced a gen- 
eralized scattering length that can be used as an alterna- 
tive parameter to characterize atomic interaction at low 
temperatures, and discussed the changes of threshold be- 
haviors in the special cases of zero and infinite generalized 
scattering lengths. The analytic formulas derived make 
it possible for an accurate description of atomic interac- 
tion in the ultracold regime without having to know the 
details of the QDT or resort to numerical calculations. 

The generalities of the results of this work will further 
manifest themselves in a companion work on atomic in- 
teraction around a magnetic Feshbach resonance, which 
is necessarily multichannel in nature P, H, Q . We will 
show that such a multichannel problem can be rigor- 
ously reduced to an effective single channel problem, to 
which most of our results here remain applicable. The 
key difference will be that the effective Ki'^ parameter 
becomes energy dependent. It is this energy dependence 
that leads to deviations from the single channel universal 



behaviors that correspond to the results here specialized 
to an energy-independent -K^f*^. 
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APPENDIX: ANALYTIC RESULTS OF Kf AND 
GENERALIZED SCATTERING LENGTHS FOR 
TWO TYPE OF MODEL POTENTIALS 



In a previous work [5l|, we have derived, for two classes 
of model potentials, the analytic results of K''{e — 0,1) 
and the number of bound states Ni for an arbitrary I. 
One class, denoted by HSTn, is of the type of a hard- 
sphere with an attractive tail: 



VHSTn(r) 



oo , < f 



(A.l) 



The other, denoted by LJn, is of the type of Lennard- 
Jones LJ{n, 2n ~ 2): 



Vij„(r) = -a„/r"-f C2„-2/r 



2n-2 



(A.2) 



which corresponds, in particular, to a LJ(6,10) potential 
for n = 6. 

For HSTrt potentials, we have shown that the K'^ pa- 
rameter at zero energy is given by (sij 



-^1^0(^0) cos(7ri/o/2) - Y^g{yQ) sin(7rt/o/2) 



J„„{yo) sin(7ri/o/2) + Y^g{yo) cos{-kvq/2) ' 

(A.3) 

where vq = (21 + l)/{n — 2), J and Y are the Bessel 
functions [13], and yo [2/(n - 2)](^„/ro)("-2)/2^ in 
which /3„ = (2/iC„/?i^)^/'-"^'^' is the length scale associ- 
ated with the Cn/r" type of potentials. 
For LJn potentials, we have derived (sH 

X£j„(0, 1) = tan(^z/o/2)[l + /i/(zo)][l - hi{zo)]-\ (A.4) 

where zo = (/3„/^2«-2)""V[2(?^ - 2)] and 

^ sin Trjzo + 1/2 - iyo/2)T{zo + 1/2 - 1^0/2) 

' ^" sin 7r(zo + 1/2 + j/o/2)r(zo + 1/2 + i/o/2) ' 

(A.5) 

Substituting these results into Eq. We obtain 



iff"(e = 0) = -J.„(yo)/>^„(2/o) 
for HSTn potentials, and 

sin(7ri'o)/i;(2:o) 



K^ie = 0) 



1 - cos(7riyo)/i;(zo) 



(A.6) 



(A.7) 



for LJn potentials. Both results are applicable for arbi- 
trary n and /. 
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Specializing to the case of n = 6, these results, com- 
bined with QDT for n = 6 [l^ [2^, [5^ , provide an accu- 
rate description of the scattering and bound state prop- 
erties for these model potentials over a wide range of 
energies around the threshold. They also give us the an- 
alytic results of the generalized scattering length for an 
arbitrary I. From Eq. (j48p . we have, for n — 6, 

ai = ai[{-iy -Y,„iyo)/JM] , (A.8) 

for HST6 potential, where vq = {21 + l)/4 and yo — 
{Pe/ro)y2, and 



sin(7ri'o)/i;(zo) 

for LJ6 potential, where zq — {f3e//3io)'^/8. The gener- 



alized effective range can be derived from these results 
using Eq. ([49|l . If one specializes Eq. IjA.Sp to / = 0, one 
recovers, for the HST6 potential, the result of Gribakin 
and Flambaum [19] for the s wave. 

The analytic results of this Appendix are useful for 
designing model potentials to investigate universal prop- 
erties not only in two-body [36l. [5ll | , bu t also in few-body 
and many-body quantum systems [sil, [sl, [s^ . They can 
also be used to check the accuracies of various numerical 
techniques and methods. Last but not the least, they 
give an explicit illustration of one of the important prop- 
erties of K'^ and related parameters, i.e., K'^{e,l), and 
therefore Ki^(e = 0) and ai, are meromorphic functions 
of Z, namely, they are analytic functions of I with only 
simple poles (25l |. 
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